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Abstract. One unexpected recent result from heavy-ion collisions is the large 
' suppression and elliptic flow of electrons from heavy flavor decay. Further 

measurements of properties of electrons from heavy flavor decay are crucial 
O ' to understanding the origin of this suppression and its implications for the 

properties of the hot matter produced at RHIC. Two particle correlations have 
been used extensively to study the propagation of hard partons through the 
produced matter in heavy-ion collisions. Measurements in p+p collisions are 
important both as reference for heavy ion measurements and to study heavy 
flavor production and fragmentation in the vacuum. Preliminary results of 
correlations between electrons from heavy flavor decay with charged hadrons 
from p+p collisions are shown. 
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1. Introduction 



$_i ' Two particle azimuthal correlations at the Relativistic Heavy Ion Collider have been 



shown to be enormously powerful to study the propagation of jet partons through 
the hot nuclear matter produced in nucleus- nucleus collisions. The same measure- 
ments in proton-proton collisions provide an essential baseline for interpreting heavy 
ion results and have also been interesting in understanding jet fragmentation [[1]. 

As of now these jet- induced correlations have only involved light hadrons. How- 
ever, with recent results showing the suppression of electrons from the decay of 
charm and bottom mesons at a level similar to that of tt° [ [2] there is great interest 
in looking at the correlations of these electrons with other hadrons in the event 
to study their jet structure and how it differs from that observed in proton-proton 
collisions and also from the jet structure of light hadron jet correlations. The proton- 
proton collisions are especially necessary for heavy flavor correlations because the 
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measured electron is a decay product, thus hadrons from the heavy meson decay 
contribute to the correlations in addition to other fragmentation products. 



2. Analysis Method 

Studies of electrons from heavy flavor decay are complicated by the large back- 
ground of electrons from light meson decay and photon conversions. In this anal- 
ysis we statistically subtract the correlations from these photonic sources using 
a method similar to that used to measure direct photon-hadron correlations [ [5] . 
What is measured is the inclusive electron-hadron conditional hadron yields per 
trigger electron, Y einc -h- This is a weighted average of conditional yields hadrons 
associated with electrons from electrons from heavy flavor decay and the conditional 
yield of hadrons associated with electrons from photonic sources: 

yr _ Ne HF Y eHF -h + N ephot Y £phot -h . . 

r e inat -h — AT , AT "J 



Ne HF + N ephot 



Y eHF -h is then given by: 



v _ {Rhf + l)y e ,„ e -fe - Y ephot - h 

le HF -h — 5 {*) 

Khf 

where Rhf = jv" HF an< ^ ^ s taken from the published PHENIX measurement [0] 
for the electron p T bins used in this analysis. 

The remaining unknown is Y e hot —h- I* 1 the electron px range of interest here, 
1.5< px <4.5GeV/c the photonic electrons are dominated by tt° Dalitz decay and 
photon conversions [ |4j , where the primary source of the photons which convert is 
7r — » 77 decay. We then measure inclusive photon-hadron correlations, where the 
inclusive photons are dominately from ir° decay. Simulations are used to construct 
Y e phot -h from Y lmc - h : 

Y e phot -h(PT,i) = W i(PTd) Y l^c-h(PT,j) (3) 



where each i and j represents a bin in electron or photon px, respectively, 0.5GeV/c 
wide. 

Two methods are used to determine the Wi(j>T,j) coefficients. In the first method 
the measured inclusive photon spectrum in p+p collisions is input into a GEANT 
based simulation of the PHENIX detector. The same electron identification cuts 
as in the real data analysis are then applied to reconstructed conversion electrons 
and the relationship between the input photon px and the reconstructed conversion 
electron px determines w. In the second method, the measured ir° spectrum from [ 
[S] is taken as input to a Monte Carlo simulation which decays the ir°s via their 
Dalitz decay. The relationship between the intermediate off-shell photon and the 
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Fig. 1. w for electrons with 3.0< px <3.5GeV/c from the method using Dalitz 
decays (black solid line) and photon conversions (blue dashed line) as a function of 
the photon px- 



resulting electron are used in the same manner as in the first method to determine 
w. The resulting weights are shown in Fig. [T] for a single electron px selection. The 
difference between the two methods is small. The ir° spectrum falls steeply with 
Pt and thus the measured photonic electrons at a given px are dominated by those 
carrying a large fraction of the ir° px regardless of whether the electron comes from 
a conversion or a Dalitz decay. 

The measured conditional yields are corrected for the hadron acceptance and 
reconstruction efficiency by comparing the raw measured hadron px spectrum to the 
corrected spectrum in Ref. The non- uniform azimuthal acceptance in PHENIX 
is corrected for by using mixed pairs. The jet functions are obtained by subtracting 
the combinatoric background by the zero yield at minimum (ZYAM) [[7J method. 

The electrons were identified as in Refs. [ HI |5]. The systematic errors are 
from the different methods used to construct Y ephot -h, the statistical error on the 
determination of the ZYAM level which moves all points by the same amount in 
a given px, e and px,h combination. For Y eHF _h there is an additional uncertainty 
from the Rhf value used to subtract the jet functions. These systematic errors are 
shown the grey boxes in the plots that follow. The systematic error on the efficiency 
correction is 10% and moves all conditional yield measurements together and is not 
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Fig. 2. Efficiency corrected ei nc — h jet functions (solid black) and e p hot — h jet 
functions (hollow red) conditional yields for electrons with 2.0< pr <3.0GeV/c and 
seven hadron pr bins from 0.5-4.5GeV/c. 



shown on the plots. 
3. Results 

Fig.[2]shows the corrected jet functions for inclusive electrons with 2.0< px <3.0GeV/c 
and seven hadron pt bins. The constructed photonic electron jet functions are 
shown in the hollow points. The fits shown are to near and away side Gaussians 
with a flat background. The resulting e#F — h jet functions are extracted using 
Eq. [3 and are shown in Fig. [3) 

In heavy quark fragmentation most of the quark momentum is carried by the 
heavy meson [|H |TT], thus the near side eHF ~ h is expected to be dominated by 
pairs in which both particles come from the heavy meson decay. Due to the decay 
kinematics, we then expect the Gaussian widths of the near side to be wider for 
p-hf — h than ei nc — h. The widths from both categories of electrons are shown in 
Fig. [U The measured near side widths for chf — h are wider than e p hot — h and 
are in agreement with widths from charm production in PYTHIA [ 14 . 

The jet functions are then integrated to get near and away side conditional 
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Fig. 3. Y eHF — h conditional yields following the subtraction in Eq. [2] for 2.0< 
PT, e <3.0GeV/c and seven hadron pt bins from 0.5-4.5GeV/c. In this pr bin 
Rhf =1.1. 



yields of hadrons per trigger heavy flavor electron. The near side integration range 
is 0< A(f) <1.25 rad and the away side integration range is 1.67< A<fr < tt rad. The 
conditional yields are shown in Fig. 

The left panel of Fig. [S] shows the away side conditional yields for the data as 
a function of PT,hadron/pT,e along with exponential fits. For comparisons we also 
show the same distributions and fits for a PYTHIA simulation of charm production 
(MSEL=4). The negative slopes are shown from the data and PYTHIA are shown 
in the right panel of Fig. [5J The ratio Pt, hadron I PT,e is defined as zt for the light 
flavor jets. We can compare the away side slopes measured here to those measured 
in Fig. 19 of Ref. [[1] for ir° — h correlations. At the same pT,trig the eHF — h slopes 
are harder than the 7r° — h. However, in both cases the PT.trig is not the PTparton- 

Since the current measurements are at low pt where sizable contributions of 
electrons from charm quarks the PYTHIA comparisons made here are only for 
charm production. In order to make more detailed comparisons, it is necessary to 
include bottom quarks as well in addition to including next to leading order effects 
not included in PYTHIA. A further complication is that the heavy flavor quark that 
leads to a trigger is not necessarily balanced by an away side balancing heavy quark. 
Ref [[T5] shows a leading order calculation of the subprocess leading to mid-rapidity 
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PHENIX PRELIMINARY, Run 5 & 6 p+p\ls = 200GeV 
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Fig. 4. Gaussian widths as a function of pr,h for four electron px bins. Circles show 
the euF — h widths and open diamonds show e p hot — h widths. Pluses show euF — h 
widths from PYTHIA and the triangles show an earlier PHENIX measurement of 
chf — h widths with lower statistics and 1.0< pt <4.0GeV/c. 



D production and finds that in the px range relevant here gg — > cc contributes «20% 
and the rest of the contribution is from processes such as eg — > eg or cq{q) — > cq(q). 
These processes are included in next to leading order calculations where there are 
three partons in the final state [1 16], However the PYTHIA simulations here only 
include gg — > cc. Further studies using NLO Monte Carlo calculations are necessary 
to understand the partonic contribution to the away side jets. 

4. Conclusions and Outlook 

Jet physics with heavy flavor is an important new frontier in heavy ion physics. 
Heavy flavor triggered two particle correlations provide complementary informa- 
tion to 7T° and direct photon triggered measurements. However the measurement 
are difficult because of the background from photon conversions and light meson 
Dalitz decays. We have shown a new method which allows the extraction of correla- 
tions triggered by electrons from heavy flavor decay and showed preliminary results 
from proton- proton collisions. These measurements are an essential baseline for 
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Fig. 5. Near (left) and away (right) side conditional yields for e#F — h correlations 
for four electron pt bins and eight hadron pt bins. 




Fig. 6. (left)Away side conditional yield measurements as a function of 
Pt, hadron / Pr,e- Fits for the data and the PYTHIA are to exponentials, (right) 
Negative slopes from the data and PYTHIA fits. 
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understanding heavy ion results. These measurements are crucial to understand- 
ing the interactions between heavy quarks and the hot nuclear matter produced in 
relativistic heavy ion collisions. 

Current measurements are limited by statistics and by the inability to separate 
the correlations from charm and bottom quarks. Silicon vertex detector upgrades 
and higher luminosity data will greatly help these measurements. Additionally, 
these types of measurements will be very interesting at the Large Hadron Collider. 

Preliminary results presented at the conference on extracting the relative con- 
tributions to the electron spectrum from c and b quarks have been superseded by 
final results, submitted for publication [[8]. 
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